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Abstract 



We present ground-state calculations for laterally coupled quantum dots con- 
taining 2, 4, and 8 electrons. As our emphasis is on spin effects our results 
are obtained by applying spin-density functional theory (SDFT). By varying 
the distance between the centers of the coupled quantum dots, the transition 
from weak to strong coupling situation is realized. For the 2-electron system 
we also apply the Heitler-London approximation and analytical concepts to 
check the reliability of SDFT calculations in this case. In addition we dis- 
cuss the features of the Coulomb staircase of laterally coupled quantum dots 
in the weak and strong coupling regimes in comparison to that of a circular 
parabolic quantum dot. 
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I. INTRODUCTION 



Semiconductor quantum dots are manmade nanoscale structures in which electrons are 
confined in all three spatial directions similar to the physical situation in atoms. As they 
show typical atomic properties like discrete energy levels and shell structures they are often 
referred to as artificial atomsB. However, in contrast to natural atoms, in quantum dots the 
number of electrons is tunable and the characteristic lengths of the system corresponding 
to external confinement potential, electron-electron interaction, and an applied magnetic 
field are of comparable size. Therefore, these systems are ideal objects to study interaction 
effects such as transitions in the ground-state spin configuration as a function of the magnetic 
field. 

Starting from quantum dots as a structure which is well understood by now, more com- 
plex systems are conceivable and likely to have perspective in future applicationsa'Ei'D. A 
simple example is analogous to a two-atom molecule consisting of two coupled quantum 
dots. The more principal properties of coupled quantum dots follow from the simple model 
of e.g. two disk-like dots side by side in a plane or on top of each other (lateral or verti- 
cal quantum-dot molecules (QDMs)), respectively, which are the subject of experimentafli'B 
and theoreticali'i'll2l'0 studies. Tunneling and overlap of orbital wavefunctions lead to bond- 
ing and antibonding states with defined parity, which are delocalized over the two-dot system 
as in natural systems. The advantage of QDMs is the tunability of the electron number and 
of the coupling strength by proper design. The latter allows to investigate the transition 
between weak and strong coupling. 

The increased interest in coupled quantum dots is indicated by a growing number of 
papers on this topic: Austing et alB measured addition spectra of vertical QDM in different 
coupling regimes. The characteristic features found in their spectra have been confirmed 
by Imamura et a/.i and Rontani et a/.i applying exact diagonalization techniques. These 
ground-state calculations were performed for up to 6 electrons as a function of the interdot 
distance and an external magnetic field. Coulomb blockade effects in lateral QDMs have been 
measured by Waugh et a/.&B, who found a conductance pattern similar to that of single dots 
for strong coupling, whereas the weakly coupled systems showed a pairing of conductance 
peaks. Theoretical work on ground-state properties of lateral QDMs was performed by 
Yannouleas and Landman0 using an unrestricted Hartree-Fock approach and by Nagaraja 
et alB applying SDFT. 

In this paper we investigate the electronic structure of the ground state of two laterally 
coupled (identical) quantum dots as a function of the interdot distance. In our calculations, 
applying SDFT, we focus on spin-dependent effects in few-electron quantum dot molecules 
containing 2, 4, and 8 electrons. In the case of 4 (8) electrons we find a change of spin 
configuration with decreasing coupling strength from = I {Sz = 1) to S'^ = {S^ = 2). 
For 2 electrons SDFT yields a transition from spin singlet to spin triplet. It turns out, 
however, that the latter result is an artifact of the SDFT, as can be shown (at least for 
our choice of the confinement potential) by using alternative approaches. Furthermore, we 
show the characteristic Coulomb staircases for the chemical potential in the weak and strong 
coupling limits, respectively. 

For all our calculationslli we assume that the two dots are coupled in a quantum- 
mechanically coherent way. This means we use delocalized electron states extended over 
the whole quantum dot molecule whereas states localized to the left or right dot do not 







exist in this description. Note, that this model does not contain all aspects of the actual 
physical situation where an electron occupying a covalent state interacts with other elec- 
trons and charges in the surrounding semiconductor material. These interactions may cause 
a dephasing of the quantum mechanical wavefunction leading to a breakdown of the delo- 
calized stated. Up to now no theory is available which allows a reliable approximation for 
dephasing rates in realistic nanostructures. 

The outline of this paper is the following: In the next section we describe our model of the 
QDM and discuss the single-particle spectrum of the non-interacting system. In Sec. Ill we 
sketch the basics of SDFT and of the numerical method for the simultaneous solution of the 
Kohn-Sham (KS) equations and calculation of the selfconsistent potentials. The following 
part is devoted to the 2-electron QDM. We compare the SDFT results with exact (analytical) 
results and with the Heitler-London approach to point out a particular problem of SDFT. 
Sec. V relates the spin configuration for 4 and 8 electrons to the single-particle spectrum and 
summarizes the corresponding SDFT results. Finally we show the typical Coulomb staircase 
of the addition energies for the weak and the strong coupling limits, respectively. 



II. THE MODEL AND ITS SINGLE-PARTICLE SYSTEM 

We consider quantum dots prepared by laterally confining two-dimensional electrons in a 
GaAs/AlGaAs heterostructure (material parameters of GaAs: m* = 0.067me, e = 12.4). As 
the lateral confinement is much weaker than that in growth direction (z-direction) we adopt 
the standard situation of the electronic quantum limit and describe the electron density 
by n{x,y,z) = n{x,y)S{z). Thus, our system is defined by the effective two-dimensional 
Hamiltonian 



(1) 




j=ifc=j+i |rj ffel 



with pj and being vectors in the xy-plane. Our choice of the external confinement 
potential V{r) is 

V{r) = ]^m*ujl min ((r - L)' , (r + L)') , (2) 

where the coupling strength depends on the distance d = |2L| between the potential minima 
at ±L in the xy-plane (Fig. 1). Along the line between these minima there is a barrier of 
height 

VM = \^*u^l(^^' ■ (3) 

This model potential describes a single parabolic quantum dot for = as the extreme case 
of strong coupling and two separate quantum dots of the same shape in the weak-coupling 
limit for d ^ oo. Due to the broken axial symmetry of V^(r) for d 7^ the numerical 
requirements are increased in comparison with those for a single circular quantum dot. 

In Fig. we show the lowest energy levels of the noninteracting electrons in dependence 
on the distance between the dot centers for the typical confinement energy Kuq = 3 meV. 



For d = 0, when the confinement potential of the QDM degenerates to the simple parabolic 
potential with axial symmetry, the spectrum shows the typical shell structure of an isotropic 
harmonic oscillator. For large distances d we obtain the same spectrum, now twofold because 
the system consists of two identical but completely separated quantum dots. By switching 
on the coupling (or lowering d starting from the weak-coupling limit) we recover the antic- 
ipated properties of a diatomic molecule: the energies decrease (increase) with respect to 
the reference level due to formation of bonding (antibonding) states, which have even (odd) 
parity as indicated by solid (dashed) lines. The intermediate coupling regime is dominated 
by crossings leading to the rearrangement of the level structure. In the weak-coupling limit 
the nth bundle of the twofold level structure with energy E ^ nhuo consists of n bonding 
states and n antibonding states. 



III. SOFT AND KS EQUATIONS 

In order to include the electron-electron interaction we employ the spin-density functional 
theory (SDFT). It is the generalization of the DFT formalism, originally established by 
Hohenberg, Kohn, and Sham0'0, to spin-polarized systems^ by including the coupling of 
the magnetization to an applied magnetic field. Accordingly, the Hohenberg-Kohn (HK) 
theorem has to be modified with respect to the spin degrees of freedom^. For this case, 
it states that two different nondegenerate ground-state wavefunctions |\E') and l^'') always 
yield different tupels (?T-(r), m(r)) 7^ {n'{r), m'(r)) of electron density n{r) and magnetization 
m(r). This is sufficient to establish a functional of the total energy with the usual functional 
properties 

Evo,Bo[n, m] = FhkK m] + y" dr [V{r)n{r) - B(r) ■ m(r)] (4) 

and the universal HK functional 

FhrKhi] = {^[n,m]\T + W\'^[n,m]). (5) 

In the limit B — > considered here, the SDFT scheme can yield a spin-polarized ground 
state for even electron numbers due to Hund's rule0'lli. These effects have already attracted 
much interest and will also be discussed in this paper. 

The spin-degree of freedom is considered in the Kohn-Sham (KS) equationsEl by assuming 
the total spin Sz in z-direction to be a good quantum number 




+V^ci[n+, n.], r] U-(r) = e^^r) (6) 



with the spin cr = ± in 2;-direction and the KS energies < £2 — ••• • ^ system 
containing particles we calculate the occupation numbers of the KS levels in the ground 
state due to 

7j = 1 < /i (7) 

< 7j < 1 ej = fi (8) 
7 • = e,>fi (9) 
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(yu: chemical potential) with the constraints 

j:ij=n^ (10) 

N+ + N~ = N. (11) 

This leads to the spin densities 

W = E7j V'Jlr)', (12) 



■Tin 



the total density 

n(r) = ri+(r) + ri_(r), (13) 
and the magnetization in z-direction (/i^: Bohr's magneton) 

^^(r) = -/iB (^+(r) - n_(r)) . (14) 
The exchange-correlation (xc) potentials 

(\r> n ] r)- ^-^xcK,n_] 
Vxc(K,n_J,r)- ^^^^^^ (15) 

are defined as functional derivatives of the xc energy functional 

1 /■,/■, ,n(r)n(r') 



i e f f 1 
Exc[n+,n_] = FuK[n+,n^] - i^-^;^ j '^^ j dr'- 



'0 

-rsK,n_]. (16) 

(Ts[n+,n_] denotes the kinetic energy functional of the KS system.) The total ground-state 
energy Eq of the interacting system can be calculated from 

i?o = E7r4--^/dr/dr'-"^^^"^^'^ 



^.^ ^ ^ 2 47re£o-/ ^ |r - r' 



+ExcK,n_] -E / drVx"c(K,^-],r)n.(r). (17) 

Concerning the xc potentials we apply the local spin density approximation (LSDA) 

Exc[n+,n_]^ y dr (n+(r) +n_(r))£xc(^^+(r),n_(r)) (18) 

and the Pade approximation of £:xc(^+(r), 'n.-(r)) in two dimensions following Tanatar and 
CeperleyEl. 

For the numerical solution of the KS equations we calculate all quantities with spatial 
dependence on a two-dimensional grid in real space. The damped gradient iteration^ ensures 
a simultaneous solution of the KS equations together with the corresponding self-consistent 
potentials. This method uses the iteration scheme 

I^C^+i)) = (19) 



with the KS Hamiltonian H (only positive eigenvalues), which converges to the ground-state 
wavefunction. The inversion of the Hamiltonian is performed approximately0'@ using 




instead of Eq. ([T9|), where xq, -E^o, Eyo are iteration parameters and , the asymptotic 
contributions of the external potential for |r| oo. Excited state KS wavefunctions are 
calculated using the same iteration scheme with an additional orthogonalization routine. 
The Hartree and xc potentials depending on the densities of the KS wavefunctions are 
recalculated in each step. Thus, the KS Hamiltonian is modified in each step of the iteration. 



IV. RESULTS FOR THE TWO-ELECTRON QDM 

Fig. 1^ depicts the KS energies for a quantum dot molecule containing two electrons as a 
function of the distance d between the centers of the molecule (we assume Uujq = 3 meV). 

Starting with d = the QDM degenerates into a circular parabolic quantum dot showing 
the known shell structure of the noninteracting case. Due to the closed shell for 2 electrons 
the ground state is not spin-polarized. This ground-state spin configuration is stable for 
increasing d as long as the energy gap between the levels (a) and (b) in the single-particle 
spectrum (Fig. |^) is not too small. At d ~ 4.5 ag (ap = 9.79 nm denotes the effective 
Bohr's radius for GaAs) the unpolarized ground state changes into a polarized one. This 
state with aligned spins persists as SDFT ground state even for d — > oo, because the energy 
gap between states (a) and (b) continues to decrease with increasing d. This result, similar 
to that obtained e.g. by Nagaraia et a/.li^, differs from those of the physically analogous 
problems of the vertical QDMsi'i^ and from the hydrogen molecule whose ground states are 
spin singlets. 

The SDFT result is also in contrast to a mathematical theorem stating that the spin 
configuration of the ground state of a two-electron quantum dot molecule, with our definition 
of the external confinement potential, should be a spin singlet state: As the total Hamiltonian 
H, Eq. (|l|), does not depend on spin coordinates we can focus on the spatial wavefunctions. 



Due to Theorem XHI.47 of Ref. ^ the (spatial) ground state wavefunction of the considered 
two-body-problem is positive and nondegenerate, this means it is symmetric in the spatial 
coordinates. As a consequence of Pauli's principle the ground state has to be a spin singlet 
state. 

In view of this rigorous statement, the SDFT results of Fig. |^ for d > 4.5 cannot 
be considered as true ground state values (the same holds for the corresponding SDFT 
densities) . 

We test this finding by using a perturbation approach analogous to the Heitler-London 
model for the hydrogen molecule, which in contrast to SDFT has the advantage of yielding 
an upper limit for the ground state energy. Similar calculations for other shapes of the 
external potential were done in Ref. 



For the two-electron QDM we rewrite the total Hamiltonian (0) as 
^ = ^ + ((ri - L)' > (ri + L)'' 

\ + o"^*^o n^in ((r2 - L)^ , (r2 + L) 



2m* 2 
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47r££:o |ri - 

Hf + /J2^ + t/f + U^ + W (21) 



with 



<'' = ^ + ^-*-o^(r.TLf (22) 

Uf^^ = ]^ni*ul mill (0, ±AYj ■ L) (23) 
1 

W = -—- -. (24) 

A-keeq |ri - r2| 

The ground-state wavefunctions of the shifted harmonic oscillators H^l^ are denoted by 

|\E'^'''"). So we can use the typical Heitler-London ansatz for the (unnormalized) singlet 
\^^) and triplet \^~) spatial wavefunctions 

|^±) = 1^ (|^R)(i)|$L)(2) ± |^R){2)|^L^(1)^ . (25) 



to calculate the expectation value of the ground-state energy of the system 



Tp± _ 



{l)(^R|(2)(^L||^|^R)(l)|^L){2) 

^ l±e-2^V/^ 
^ l±e-2iV/^ 

(l)(^R|(2)(^L|y^|^L)(l)|^R)(2) 

± 1 ± e-2^V/^ • (26) 

Note that all matrix elements are real quantities. Ih = denotes the characteristic 

oscillator length. The result for the matrix elements is 

= 2huJo f--^e-^'/'^ + ^ (1 - erf {L/h))\ (27) 

(1)(^R|(2)(^L|P^|^R)(1)|^L)(2) = 

= Uu,J^K-^^l^ll,{L-lll) (28) 

V L cIq 
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V 2 cLq 

with the modified Bessel's function lo{x) and the error function erf (x)ii. 

In Fig. ^ we compare the ground-state energies resulting from SDFT calculations with 
the Heitler-London approach. For < d < 4.5 a.Q the SDFT singlet ground-state energies are 
lower than the HL energies. As the comparison between SDFT and exact diagonalization 
for 2-electron quantum dots shows good agreementil, we expect this to hold also for the 
QDM with small d. Moreover, the HL triplet energy diverges for d ^ and becomes larger 
than the HL singlet energy. Therefore, we regard the SDFT results for small ci as a better 
approximation of the ground state than the HL results^. 

If, however, d is increased beyond 4.5 ag the HL energies become lower than the SDFT 
results and replace them as ground-state energies. It can be shown analytically that for 
sufficiently large d the HL energies fulfill the relation E^^i^ < -Ehl- Thus, in the limit d oo 
the HL results are consistent with the mathematical theorem quoted above. In addition, in 
the weak-coupling limit the interaction is reduced to the classical Coulomb repulsion between 
point charges and the total energy approaches asymptotically E = 2hu!Q + e"^ / {Aireeod) . This 
is also consistent with the energies of the HL approach. 

The energy difference of 0.36 meV between -E^^px -^hl weak-coupling limit 

can be ascribed to the self-interaction within the SDFT scheme: The charge density of 
a weakly coupled quantum dot molecule is mainly localized in its centers (charge 1 e~ per 
center) . However, within the SDFT formalism the energy of each charge in its own Coulomb 
field is counted additionally leading to an enhanced total energyil. To illustrate this, let us 
to apply the SDFT concept to a single-particle problem. For a circular parabolic quantum 
dot with hujQ = 3 meV this gives a ground-state energy of 3.18 meV, i.e. the self-interaction 
error is 0.18 meV. This value has to be doubled for the d ^ oo limit of the QDM. Thus, we 
identify the difference E^j^-pr^ — the self-interaction energy in this limit. 

The energy difference between the singlet state and the triplet state in the SDFT ap- 
proach can be qualitatively explained as an xc effect: Both spin configurations show similar 
total electron densities. Consequently the contributions from the Hartree energy and the 
spin-independent xc energy can be taken to be approximately equal for singlet and triplet 
states. However, the triplet total energy is in addition lowered by spin-dependent xc effects 
leading to the energy difference in Fig. 

To summarize this section: for small d {0 < d < 4.5 slq) we find the SDFT (singlet) 
energies to be the better approximation to the ground state than the HL results, while for 
larger d {d > 4.5 ag) the SDFT energies become too large due to the self-interaction error 
and the HL results provide the better approximation of the ground-state energy. 

V. RESULTS FOR QDMS CONTAINING 4 AND 8 ELECTRONS 

In this section we present our results for the ground-state energy and spin-configuration 
for artificial molecules containing 4 and 8 electrons as a function of the interdot distance 



s 



d. Our calculations show that a spin-polarized ground state is accompanied by a (quasi- 
) degeneracy at the Fermi level in the noninteracting single-particle spectrum whereas un- 
polarized ground states require a sufficiently large gap between the highest occupied and 
lowest unoccupied single-particle levels of the noninteracting spectrum. 

Single dots with 4 electrons have a spin-polarized ground state {Sz = 1) due to Hund's 
rule0'lll. This is the case also for the QDM with d = (Fig. 5). The spin-polarized ground 
state persists with increasing d as long as the xc energy overcompensates the energy costs of 
putting electrons with parallel spins into the bonding and antibonding states ((b) and (c) in 
Fig. ^ evolving from the second level at d = 0. For d > 2 it becomes energetically more 
favorable to occupy both the bonding level (a) and antibonding level (b) (Fig. ^ with two 
particles forming an 5*^ = ground state. This configuration remains stable even in the weak 
coupling limit and reflects the noninteracting particle picture (Fig. ^ which does not show 
level crossings at the Fermi energy for = 4 and large d. This result is in agreement with 
the work of Nagaraja et a/.ini (laterally coupled dots), Rontani et a/.i, and Imamura et a/.i 
(vertically coupled dots). On the other hand, Nagaraja et al^ do not find a spin-polarized 
ground state in the strong coupling limit, because their confinement potential obtained by 
simulating external gates does not approach that of a circular parabolic dot as in our case. 

Although even in the limit of large d a quantum-mechanical coupling (delocalized bonding 
and antibonding states) between the dots can be observed within the SDFT scheme the 
resulting density (Fig. ||) and energy (Fig. |]) tend towards a semiclassical picture: the 
densities of the dots are well separated and the total energy of the two coupled dots consists 
of twice the (quantum-mechanical) energy of a single dot containing 2 electrons and the 
energy of classical repulsion of two point-like charges (2e~) in the distance d (see inset of 
Fig. H). In contrast to the 2-electron QDM, this test for the asymptotics of the total energy 
being fulffiled underlines the reliabilty of the SDFT results for 4 electrons. 

Turning to QDMs with 8 electrons (Fig. 0) we find for d = the expected spin-polarized 
ground state {S^ = l)0'lli. Increasing d up to 7.5 this state persists with double-occupancy 
of the levels (a), (b), and (c) of Fig. ^ and single-occupancy of the quasi-degenerate levels 
(d) and (e) with aligned spins. At d ^ 7.5 ag the levels (c), (d), (e), and (f) of Fig. 
get close enough that an alignment of all four spins at the Fermi level can lower the 
total energy. Therefore, we find a spin-polarized ground state in the weak coupling regime 
(d > 8 ag) in agreement with Nagaraja et a/.0. Although the spin-aligned ground state 
was ruled out in the case of 2 electrons we believe that the Sz = 2 ground state for the 
8-electron artificial molecule in the weak-coupling regime is reliable. This is emphasized by 
the correct asyniptotics of the total energy in the limit d ^ oo (see inset of Fig. |^. Exact 
diagonalizationsi"! for vertically coupled dots performed for up to 6 electrons yield a spin- 
polarized ground state {S^ = 1) for the 6-electron artificial molecule in the weak coupling 
regime, too, indicating that spin-polarized ground states can exist in principle for QDMs. 

Also for N = 8, in the weak-coupling limit, the total energy and density can be under- 
stood in a semiclassical picture: the densities (for sufficiently large d) are well separated 
(Fig. H) and the total energy follows that of two identical dots each containing 4 electrons 
plus the Coulomb repulsion between two point-like charges (4e~) in distance d. However, 
besides the quantum mechanical phenomenon of delocalized bonding and antibonding states 
we find a (nonclassical) xc energy related effect leading to the spin-polarized ground state. 

Yannouleas and Landman^ find from unrestricted Hartree-Fock calculations electron 
localization in the form of electron puddles and Wigner supermolecules for a ratio A between 
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the oscillator length and Bohr's radius of about 1.4. In contrast, (exact) Quantum Monte 
Carlo calculationsEa indicate a transition from the Fermi liquid to the Wigner crystal regime 
if A is increased beyond a critical value of 4. Thus, our calculations performed for A ~ 2 are 
obviously in the Fermi liquid regime. 



VI. CHEMICAL POTENTIAL 

Although recent experiments succeeded in identifying molecular properties of coupled 
quantum dots like bonding and antibonding stateslll the most convenient comparison be- 
tween theoretical and experimental results is based on the chemical potential of the system. 
The chemical potential is defined as difference of total energies fi{N) = E{N) — E{N — 1) 
which can be measured by capacitance spectroscopyi'&il. 

As it turns out, the dependence of N on fi (Coulomb staircase) exhibits features which 
depend on the coupling strengthi'i'0. While in the strong coupling limit the capacitance 
spectrum of the QDM is similar to that of the single quantum dot, it exhibits small spacings 
between the chemical potentials for = 1, 2, iV = 3, 4, iV = 5, 6 etc. for the weakly coupled 
QDM which correspond to the pairing of peaks observed in capacitance spectroscopyi'S. 
Explanations for this phenomenon provided in the literaturei are based on states localized 
at the centers of the molecule and neglect quantum mechanical coupling. Using these as- 
sumptions we can easily estimate the dependence of fi on N for our system: If the QDM 
contains an even number N of electrons one may assume that N/2 electrons are localized at 
each dot. The total energy of the artificial molecule Eu{N) can be estimated by the energy 
of the single dots containing N/2 electrons {Ed{N/2)) and the Coulomb repulsion between 
the dots which is approximated as Coulomb energy between point-like charges 

/A^\ e'^(N/2)'^ 

E„iN) = 2E.{j)+-^,Ne.en. (31) 

For odd one assumes {N + l)/2 electrons to be in one dot and (A^ — l)/2 in the other one 
and the approximation for the total energy yields 



+fm2M^'^„M. (32) 

Using these equations we can calculate the chemical potential fJ,M{N) = Eu{M) — Eu{N — 1) 
of the QDM for even and odd A^ 

MiV) = Bo(f)-^;o(f-l)+|^.iVeven. (33) 
MA')^^om-Bor^) + 5^<^^.A'odd. (34) 
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Consequently, the addition energies which are necessary to add a further electron are 

in 



A/iM(iV) = (^^^) - /^D (^^^) ' ^ > 3 odd. (36) 

These equations indicate that in a weakly coupled quantum dot with strongly localized 
wavefunctions the addition of an electron to a configuration with an even electron number 
is energetically more expensive as it corresponds to an addition of an electron to a single 
quantum dot. On the other hand, adding an electron to an odd number of particles is 
relatively cheap because the addition energy only depends on the interdot Coulomb repulsion 
which is relatively small for large d and gives rise to the typical pairing of conductance peaks. 
With growing electron numbers this model fails as the increasing Coulomb repulsion caused 
by more extended charge distributions is not included. This effect destroys the splitting of 
the paired peaks. 

A conceptional difficulty of this model is that it is based on localized states yielding 
an asymmetric charge distribution for odd electron numbers. In our SDFT calculations, 
however, the delocalized states also lead to the pairing of the conductance peaks in the weak- 
coupling regime (Fig. |^) while providing symmetric electron densities for all A^. Therefore, 
we would like to emphasize, that the pairing of peaks does not depend on electron states 
localized at one center of the QDM, but can be caused by delocalized wavefunctions as well. 
In addition, delocahzed states provide the advantage that it is not necessary to determine a 
transition from a coherent to an incoherent regime when the distance d between the centers 
is continuously increased. 

VII. CONCLUSIONS 

We have studied the ground-state properties of QDMs containing 2, 4, and 8 electrons as 
a function of the distance d between the centers of the molecule using SDFT. This concept 
includes spin effects in addition to the classical Coulomb repulsion. The results obtained 
for 4 and 8 electrons are reliable over the whole range of d and can be predicted from the 
single-particle spectrum of the noninteracting system for not too weak confinement. For 4 
electrons we find a transition from a spin-polarized to a spin-unpolarized configuration with 
decreasing coupling strength. The ground state of the 8-electron QDM can be characterized 
by a spin-polarized ground state (2 parallel spins) for strong coupling and a spin-polarized 
configuration with 4 aligned spins in the weak-coupling limit. In contrast, the two-electron 
QDM is well described by SDFT only in the strong-coupling regime whereas for large d a HL 
approach yields better results. Finally we have discussed the Coulomb staircase diagram 
for the weak and the strong coupling regime and identified their characteristic features 
considering the background of delocalized states in an artificial molecule. 

This work was funded by Deutsche Forschungsgemeinschaft (Grants No. SFB 348 and 
Ro 522/16). 
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FIGURES 



FIG. 1. Model potential for a QDM. In this picture the interdot distance is d = 13 ag and the 
strength of the parabolic potential near the centers is huo = 3 meV. 

FIG. 2. Lowest levels in the spectrum of the noninteracting system as a function of the interdot 
distance d. The levels are classified as bonding and antibonding states with even and odd parity, 
respectively. 

FIG. 3. Lowest KS levels for the 2-electron QDM as a function of the interdot distance d (the 2 
lowest levels are occupied). In the shaded area at d « 4.5 ag SDFT calculations yield a transition 
from a spin-unpolarized to a spin-polarized ground state. 

FIG. 4. Comparison between the total energies of the SDFT calculation (-E^^prp) and of the 
Heitler-London approach (-Eg^) for the 2-electron QDM. -|- (— ) denotes singlet (triplet) energies. 

FIG. 5. Lowest KS levels for the 4-electron QDM as a function of the interdot distance d (the 
4 lowest levels are occupied). In the shaded area at d ~ 2 ag SDFT calculations yield a transition 
from a spin-polarized (Sz = 1) to a spin-unpolarized (Sz = 0) ground state. Inset: Ground-state 
energy of the 4-electron QDM (crosses). The solid line describes the asymptotics for large d in a 
quasi-classical picture as a combination of the ground state energy of two single dots containing 2 
electrons {E^^^ = 11.418 meV) and the Coulomb repulsion between two point- like charges (2 e~) 
in distance d. 

FIG. 6. Ground-state density of a 4-electron QDM for d = 10 ag. The character of two 
uncoupled dots is emphasized by the localization of the density in the centers of the molecule. 

FIG. 7. Lowest KS levels for the 8-electron QDM as a function of the interdot distance d 
(the 8 lowest levels are occupied). In the shaded area at d 7.5 ag SDFT calculations yield a 
transition from a spin-polarized {Sz = 1) to a higher spin-polarized {S^ = 2) ground state. Inset: 

Ground-state energy of the 8-electron QDM (crosses) . The solid line describes the asymptotics for 
large d in a quasi-classical picture as a combination of the ground state energy of two single dots 
containing 4 electrons (Ej^^^. = 41.064 meV) and the Coulomb repulsion between two point-like 
charges (4 e~) in distance d. 

FIG. 8. Ground-state density of a 8-electron QDM for d = 10 ag. The character of two 
uncoupled dots is emphasized by the localization of the density to the centers of the molecule. 

FIG. 9. Coulomb staircase (derived from SDFT results) for the weak-coupling regime (a) 
{d = 10 ag), for the strong coupling regime (c) (d = 3 ag), and for the circular parabolic quantum 
dot (c) (d = ag). The shaded areas in picture (a) indicate the pairing of conductance peaks. 
For higher electron numbers the pairing splits as a consequence of increasing Coulomb repulsion 
caused by more extended charge distributions. 
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